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The properties of light leptoquarks predicted in the context of a simple grand unified theory and 
their observability at the LHC are investigated. The SU (5) symmetry of the theory implies that 
the leptoquark couplings to matter are related to the neutrino mass matrix. We study the resulting 
connection between neutrino masses and mixing parameters and the leptoquark decays, and show that 
different light neutrino hierarchies imply distinctive leptoquark decay signatures. We also discuss 
low-energy constraints implied by searches for charged lepton flavour violation, studies of meson 
decays, and electroweak precision data. We perform a detailed parton-level study of the leptoquark 
signals and the Standard Model backgrounds at the LHC. With the clean final states containing a di- 
lepton plus two jets, the QCD production of the leptoquark pair can be observed for a leptoquark mass 
of one TeV and beyond. By examining the lepton flavor structure of the observed events, one could 
further test the model predictions related to the neutrino mass spectrum. In particular, 6-flavor tagging 
will be useful in distinguishing the neutrino mass pattern and possibly probing an unknown Majorana 
phase in the Inverted Hierarchy or the Quasi-Degenerate scenario. Electroweak associated production 
of the leptoquark doublet can also be useful in identifying the quantum numbers of the leptoquarks 
and distinguishing between the neutrino mass spectra, even though the corresponding event rates are 
smaller than for QCD production. We find that with only the clean channel of p + ^y+jets, one 
could expect an observable signal for a leptoquark masses of about 600 GeV or higher. 



I. INTRODUCTION 



The possibility of explaining all fundamental interactions in nature in a unified framework is one of 
the main motivations for considering physics beyond the Standard Model (SM) of particle physics. In this 
respect, Grand Unified Theories (GUTs) are one of the most appealing SM extensions, as they may allow 
one to understand the origin of the SM interactions and predict both the quantization of the electric charge 
and the weak mixing angle [1]. In these theories the unification of gauge couplings is realized at a scale 
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Mqut ~ 10 14-16 GeV. The unification of the matter fields implies the decay of the lightest baryon, the 
proton | and the generic existence of scalar and vector leptoquarks (denoted by LQ unless specified 
otherwise) [3|] that carry both lepton- and baryon-numbers, since the quarks and leptons live in the same 
enlarged representation of the gauge group. 

Because of the important role played by LQs in GUTS and other SM extensions as well as their 
unique phenomenological features, the general properties and phenomenology of LQs have been studied 
extensively [4j,|5|]. Recently, a simple GUT based on the SU(5) gauge symmetry that contains light LQs has 
been proposed in Ref. |Q]. In this theory, the Higgs sector is composed of 5h, 15h and 24h; neutrino 
masses are generated through the Type II seesaw mechanism [7]; and there exist light scalar leptoquarks 
consistent with the constraints from the gauge coupling unification and experimental lower bounds on the 
proton lifetime. A distinctive feature of this framework is that the LQ-matter Yukawa interactions are gov- 
erned by the neutrino mass matrix. Consequently, different scenarios for the light neutrino mass hierarchy 
lead to broadly distinguishable patterns for LQ interactions with matter. In this paper, we explore this phe- 
nomenology in detail, concentrating in particular on its implications for studies at the Large Hadron Collider 
(LHC). 

While several phenomenological and cosmological aspects of this proposal were studied in Refs. ||8|] 
and a comprehensive analysis of the corresponding LQ collider phenomenology has not appeared pre- 
viously in the literature. Here, we report the results of such an analysis. In doing so, we consider the 
constraints on the model implied by the results of neutrino oscillation experiments, searches for charged 
lepton flavor violation, and studies of SM-suppressed meson decays. We find that although the constraints 
from low-energy flavor physics are severe, they do not preclude the possibility of observing a statistically 
significant LQ signal at the LHC if the LQ mass is sufficiently heavy. The constraints from electroweak 
precision observables (EWPO) have a marginal impact at best on the LHC discovery potential of the model. 

In searching for the LQ at the LHC, we identify the dominant production mechanisms and the exper- 
imentally fully reconstructable clean final states. We find that the "golden signals" for LQ discovery are 
the QCD pair production of the charge +2/3 LQ (denoted <I>i), with its decays to a charged lepton plus a 
down-type quark (di, i = 1 — 3). After applying a series of judicious cuts to suppress SM backgrounds that 
would produce the same final state, we find that one could expect on the order of 1000 events with 100 fb _1 
integrated luminosity for a 400—600 GeV LQ. By examining the lepton flavor structure of the observed 
events, one could further test the model prediction for the neutrino mass spectrum. In particular, 6-flavor 
tagging will be useful in helping to distinguish the neutrino mass pattern and to probe an unknown Majo- 
rana phase in the "Inverted Hierarchy" or the "Quasi-Degenerate" scenario. Assuming that one will have 
knowledge of the neutrino mass spectrum in the near future, obtained from either searches for neutrinoless 
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double /3-decay or future long baseline neutrino oscillation studies, the possible collider discovery analyzed 
here would not only provide a crucial consistency check for the neutrino mass pattern, but also indicate a 
neutrino mass generation mechanism at a fundamental level. 

Our study is organized in the remainder of the paper as follows: In Section II we present the properties 
of the LQs in the model. In Section III the constraints on the physical couplings and parameters coming 
from neutrino oscillation experiments, rare processes and EWPO are investigated. The predictions for 
leptoquarks decays are discussed in Section IV. Taking into account the effects of neutrino mass and mixing 
we show the unique predictions for the branching fractions of all leptonic decays — > die^ and iff 2 — ► d^P. 
The impact of the Majorana phases is also studied. The pair and associated production mechanisms, cuts 
necessary for SM background suppression, and the signal observability at the LHC are investigated in 
Section V. We summarize our results in Section VI. 



II. LIGHT LEPTOQUARKS AND GRAND UNIFICATION 



A general classification of the properties of leptoquarks can be found in Ref. M\. Here, we concentrate on 
the specific model realization of Ref. |60 that consists of a simple realistic extension of the Georgi-Glashow 
model niOll in which the light neutrino masses are generated through the Type II seesaw mechanism. The 
Higgs sector of this model is composed of 5h, 24h and 15h, and the matter fields live in the 5 = {d , I) l 
and 10 = (it ,Q,e C )L representations, where l T = (v,e) and Q T = (u,d). As emphasized in Section 
H one of the main features of this GUT is the possibility of having light leptoquarks consistent with all 



constraints coming from the proton lifetime and the unification of gauge couplings (g, 
on the properties of these light scalar leptoquarks 



9D. Here, we focus 



^ = (*i,* 3 ) ~ (3,2,1/6). 



(1) 



The electric charges of the leptoquarks are Q(^i) = I3 + Y = 2/3 and Q(^2) = —1/3, respectively. The 
LQ lives in the 15h representation together with A ~ (1, 3, 1), the scalar triplet field responsible for 
implementing the Type II seesaw mechanism. Using the interactions in the model [6] one finds 



Y v 5 5 15 H D Y u 



llCia 2 Al L + V2d R l?$?e a 



(2) 



where we have suppressed the generation indices on the lepton and quark fields and Yukawa coupling matrix 
It is important to notice that the SU(5) symmetry of the theory implies that the coupling of the lepto- 
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quark to matter is defined by the mass matrix for neutrinos 

V2 va 

where va/V^ is the vacuum expectation value of the neutral component of the field A. As shown in 
Ref. (60 the leptoquark & b has negative (positive) contribution to b\ — b 2 (b 2 — b^), where hi stands for 
the different beta function coefficients for the three SM gauge groups. Therefore, the presence of this field 
helps to achieve unification without supersymmetry. If the LQ mass lies in the range of 100 GeV— 1 TeV, 
one obtains gauge coupling unification in agreement with experimental observations at the low energies. 
The Lagrangian relevant for our study is given by: 

C = (D^ b )\D^ b ) +C Y -V, (4) 

where 

C Y = - [V2d R Y v l a L e aP + h.c.) (5) 

V = Ml b ^ h + X x {^ b ){H^H) + X 2 (^ b H)(H^ b ) +A 3 ($ t 6 $ 6 ) 2 , (6) 

where H T = (H + ,H°) is the SM Higgs doublet. The tree-level masses of the leptoquarks, after the 
electroweak symmetry breaking (H°) = v /y/2, are given by: 

Ml 1 = M% b + Mvl/2 and M| 2 = + A 2 ^/2. (7) 

In principle, the dimensionless couplings \\, \ 2 can be either positive, negative, or zero 1 . Thus, we have 
no a priori prediction for which of the two LQs is the lightest 2 . In our study of the collider phenomenology, 
we will focus on the case in which Mip 2 — > 0, corresponding to A2 > 0. 

Working in the basis of physical fermions, the new LQ-matter Yukawa interactions read: 

d R Ti e L *i and d R T 2 v L * 2 , (8) 

where 

r x = V2 D ] R Y V E L = T 2 VI MNS K* and T 2 = V2 D R Y V N L = V2 B K* z V£ MNS . (9) 

Here, Vpmns is the leptonic (neutrino) mixing matrix; D R , El, and Nl transform the right4ianded down 
quark, left landed electron, and left4ianded neutrino mass eigenstates into the corresponding flavor states; 



For negative Ai and/or A2, boundedness of the potential restricts the magnitude of coupling relative to that of A3 and the SM 
Higgs quartic coupling. 

2 We note that SM radiative corrections will generate a mass splitting JTlTl — M* 2 w 106 MeV in the MS scheme. 
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and K3 is a diagonal matrix containing the unknown phases 

^ 



( e ««i 



e ia2 




D ia 3 



(10) 



The other mixing matrix, B, appearing in Eq. (© involves the product of RH quark and LH charged lepton 



rotation matrices: is B = D R E* L . It has been pointed out (9, 
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IJ] that in order to satisfy the constraints 



from proton decay, the unitary matrix B must be of the form 

/ n p*A \ 




B 



(11) 



/ 



e 

e i/32 



Once we impose the constraints from the lower limit on the proton lifetime, the LQ-matter couplings are 
dictated by the neutrino masses and mixing, as presented in detail in Appendix A. It is useful to notice 
that the Yukawa couplings scale parametrically as oc m v /v^. We summarize the interactions for the 
leptoquarks in Table. HI 



III. CONSTRAINTS ON THE PHYSICAL PARAMETERS 

In this section we discuss the constraints from neutrino oscillation experiments, rare decays and collider 
experiments on the Yukawa couplings and the physical parameters , Mq, 2 and va in this theory. 



A. Constraints From Neutrino Oscillation Experiments 

The physical Yukawa couplings of leptoquarks for the leptonic decays are given by Eq. (O. In order to 
understand the constraints from neutrino physics we start from the correlation between the neutrino masses 
and mixing angles. The leptonic mixing matrix is given by 



VpMNS 



C12C13 



C13S12 



e~ i5 s 13 \ 



-Cl2Sl 3 S 2 3e 



i<5 



AS 



C23S12 C12C23 - e'"Si 2 Si 3 S23 Ci 3 S 23 



diag( e ^/ 2 ,l,e i$2 / 2 ),(12) 



\ S12S23 - e 1<5 Ci 2 C 2 3Sl3 -C23Sl2Sl3e 1<5 - Ci 2 S23 C13C23 J 



where s 



'j 



sm l 



ij, fyj 



cos Oij, < 8ij < tt/2 and < 5 < 2ir. The phase 5 is the Dirac CP phase, and 



<I>i are the Majorana phases. The experimental constraints on the neutrino masses and mixing parameters, 



6 



Fields 


Vertices 


Couplings 


*i(i = l,2) 
Yukawa 


h = e, h = v 


Pi = (1 - 75)/2 


QCD 




igs^ipi -P2)n 
i y s 2 2 


3-point EW 




iQie(pi -P2) M 
Qi =2/3, Q 2 = -1/3 






cos 9 w K l (p 1 -P2)n 
A 'i = 4 cos^u-' x 2 2 ^-5^r 




*i(pi)*3(pa)W- 




4-point EW 




G Z Z11 - l 8(s ? + ^)^™ - 18(sHg .) 












iGzAii9liv 






iGwWii9nv 

Gwwn = Gww22 = 



TABLE I: Feynman rules for the leptoquarks Yukawa and gauge interactions. The momenta are all assumed to be 
incoming. 

at 2a level are UM 



7.3 x 10" 


" 5 eV 2 


< 


Ami 


< 


8.1 x 10~ 5 eV 2 , 


(13) 


2.1 x 10" 


" 3 eV 2 


< 


lAmlil 


< 


2.7 x 10~ 3 eV 2 , 


(14) 




0.28 


< 


sin 2 6»i2 


< 


0.37, 


(15) 




0.38 


< 


sin 2 6>23 


< 


0.63, 


(16) 








sin 2 6>i3 


< 


0.033, 


(17) 



and J2i m i < 1-2 eV. Following the convention, we denote the case Am^ = m| — mf > as the normal 
hierarchy (NH) and otherwise the inverted hierarchy (IH). 

In this section we neglect all the phases. The observed neutrino mass spectra indicate that the neutrino 
mass matrix 

Mu = VhiN S ™f a9 vkiNS> ( 18 ) 
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FIG. 1: Constraints on the leptoquark couplings Yf 1 = Tj 1 x va versus the lowest neutrino mass for NH (left) and 
IH (right) when all the phases vanish. Due to the symmetry as in Eq. 1201 . the equal couplings are also indicated in 
the parenthesis. 



presents the following patterns 



M}} 

MV > M~?,M, 



33 for NH, 
for IH, 



(19) 



and M™>M}*M}? for both NH, IH. 



where the last relation for the off-diagonal elements is due to the large atmospheric mixing angle #23- See 



for example Ref. B15I1 for more detailed discussions. 

In order to see the relationship between neutrino mass pattern and the Yukawa couplings of the lepto- 
quarks, we first introduce a dimensionful coupling matrix relevant to Yukawa coupling 

/ 1 \ 



Y 1 = T 1 x v A 



V 



M u , or Yl j = M^~ i) ' i . 



(20) 



/ 



1 

1 

Here we have neglected the phases in the B matrix. Note that the neutrino mass matrix M u is symmetric 
but Y\ is not, with the first index i for the down-type quarks and the second index j for the charged leptons. 

We are thus able to determine all the Yukawa couplings that are consistent with a given set of neutrino 
masses and oscillation parameters. To illustrate this relationship, we perform a numerical scan over the 
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FIG. 2: Constraints on the leptoquark couplings y/ 2 = T'f x va versus the lowest neutrino mass for NH (left) and 
IH (right) when all the phases vanish. Due to the symmetry as in Eq. d20l i. the equal couplings are also indicated in 
the parenthesis. 

neutrino masses and mixing angles, using a uniform, random distribution while taking into account the 
above experimental constraints and neglecting all phases. We show the scatter plots for the allowed values 
for the couplings of each lepton flavor in Figs.[U|2]and[3l versus the lightest neutrino mass in each spectrum, 
the normal hierarchy (left panels) and the inverted hierarchy (right panels). We see two distinctive regions 
in terms of the lightest neutrino mass. In the case m 1 ( 3 ) < 1CT 1 eV, 

Fig. fflfor e : Y 31 ~ Y/ 1 , if 1 for NH, 
Yi 31 >F 1 11 ,Y 1 21 forffl, 
Fig. |2]for fi : Y? 2 > Y x 12 > Y 32 for both NH and IH, 
Fig. [3]for r : Y/ 3 > Y 23 > Y 33 for both NH and IH. 

On the other hand, for rnug\ > 10 _1 eV, we have the quasi-degenerate spectrum M^ 1 « M 22 m M 33 , and 
the leading couplings in each lepton flavor are 

Y x 31 « Yl 2 ~ Yl 3 . (21) 



As for the ^ decays, we sum over the final state neutrinos since they are experimentally unobservable. 
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FIG. 3: Constraints on the leptoquark couplings Yj* 3 = Vf x va versus the lowest neutrino mass for NH (left) and 
IH (right) when all the phases vanish. Due to the symmetry as in Eq. d20l i. the equal couplings are also indicated in 
the parentheses. 

Thus the relevant couplings are written as 

3 

y 2 * = ^ K j \ 2 X v\ (i = 1, 2, 3 for down quarks d, s, b). (22) 

The allowed values are shown in Fig. [4] In this case, Kf < K>\ F 2 2 in the NH and if > Y£,Y% in the IH 
when the lightest neutrino mass is smaller than 10~ 2 eV. The explicit expressions of Y\ and Y2 obtained via 
Ti and T 2 are collected in Appendix A. 



B. Constraints from Rare Decays 



Here we discuss the constraints coming from meson decays, meson-antimeson mixing and lepton flavor 
violating processes. The most important constraints for the leptoquark parameters are from the decays 



to dileptons 



;i Jl6 1. 
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FIG. 4: Constraints on = l-Hf I 2 x V A versus the lowest neutrino mass for NH (left) and IH (right) when all 

the phases vanish. 



Mesons 


Final states 


Branching fraction 


Kl 


fi+fj,- 
e + e~ 


(6.84 ±0.11) x 1CT 9 
(0.087±8;g|I) x 1(T 10 
< 4.7 x 1(T 12 



TABLE II: Experimental Constraints on K® pure leptonic decays 



1. Kl decays 

The meson leptonic decay rate is given by 

h. i \ynr\2 \-pji\2 / f \ 2 

r ( M (f >,», «•) = ^ _±_E^E^ (EiEr _ (23) 

where 

= y^L^, k = ^A^(i, (24) 
V • 2 

For Kl decay, we take m K = 497.648 MeV, = 160 MeV and t Kl = 5.116 x 10" 8 s [16]. We list 
the most precise results for the Kl decay branching fractions B(Kl — ► X) in Table [TT1 In the case of the 
and e + e~ final states, SM expectations involve the sum of short-distance (SD) and long-distance 
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(LD) contributions The latter are dominant, as the predicted short-distance contribution gives 

roughly one sixth of the measured branching ratio. Using chiral perturbation theory to perform a model- 
independent analysis of the long-distance contributions, the authors of Ref. J19I1 find that they are nearly 
saturated by the lowest order (l.o.) term in the chiral expansion that is fixed from the measured Kl — ► 77 
branching fraction. The sub-leading long-distance contribution depends on an a priori unknown scale- 
dependent low-energy constant x(aO whose value is obtained from analysis of Kl — > rj — > /U + /i~, 
and 7T° — ► e + e~ decays 3 . Of these, B(Kl — ► M + aO i s known most precisely. The difference between its 
experimental value and the sum of the short-distance and leading order long-distance contributions 



determines the value of \ implied by this decay. 

Any non-SM contribution to Kl — > fJ> + [i~ would enter this difference and, thus, affect the extracted 
value of x- So long as the non-SM contribution is smaller than the uncertainty in Ajf(/i//), it could be 
absorbed by a shift in \ within its error bars without appreciably affecting the analysis of Kl, vt°, and 77 
decays into £ + £~. The uncertainty 5Ak is dominated by the theoretical uncertainty in B(Kl — > : 



To be conservative, we thus require that the LQ contribution to B(Kl — > be smaller than the 

magnitude of 5Ak{^Ij) given in Eq. (l26l ). 

Once the constant \ has been obtained from the other processes discussed above, it is possible to make 
a model-independent prediction for the Kl — ► e + e~ branching fraction. In this case, however, the compar- 
ison of the experimental result with the Standard Model prediction, <5A^-(ee), is dominated by the experi- 
mental error in B(Kl — ► e + e~). Consequently, we require that the LQ contribution to this channel be less 
than the average of the upper and lower errors quoted in Table JT] 

In Figs. [5] and [6j we show the branching fractions for K\ — ► e + e~, fi + pT for the NH and IH, 
respectively, taking into account the leptoquark ^1 contribution versus the lightest neutrino mass for 
M^jUa >1200 GeV- eV (left panel) and 1600 GeV- eV (right panel). Note that the decay rate is pro- 
portional to (M^ 1 «a) since each T M ~ m v /v&. We see that there exist parameter space points for 
which the leptoquark contribution is smaller than uncertainty in Ak (££) as indicated by the horizontal lines 
for M^va >1600 GeV- eV. 

In Fig. |7] we show the predictions for LQ contributions to the decay channel e~fi + with AI^^a >800 
GeV- eV. In this case, the SM prediction is highly suppressed since flavor violation arises solely through the 



A K (nfi) = B(K L - /i + /O cxp " [B(K L - ^V) SD + B(K L /u + /i-)f:°] 



(25) 



6A K (w) ^ (±0.6) x 10 



-9 



(26) 



3 The scale /i is typically chose to be the p-meson mass. 
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neutrino sector, so any observable effect would have to arise from new physics. The experimental BR limits, 
indicated by the horizontal lines in the left and right panels of Fig. |7]thus translate into direct constraints on 
the LQ parameters. We find the BR results for this channel more spread out and, thus, less constraining than 
are those for the lepton flavor conserving decays. This observation leads to an important point: Although 
the experimental bounds on the e/i and ee channels are stronger than that of they actually provide less 
stringent constraints on the model parameters, due to the suppressed couplings governed by the neutrino 
oscillation data. 

Additional constraints follow from other meson properties and low-energy semileptonic interactions, 
such as the mass difference between Kl and Ks- Amx = 3.48 x 10 -15 GeV Jl^l . The box diagrams for 



a leptoquark with couplings to leptons give a contribution to ArriK 1171 



1811 . that yields a constraint 



(rir) 21 ) < 7^ -^5.2xlO- 10 M| r (27) 

' lK m K 

As we have shown in the previous section when the lightest neutrino mass is smaller than 10~ 2 eV, (Y\Y^) 21 
is smaller than 25(24) x 10 -4 eV 2 in the NH (IH) spectrum. Using these results one finds 

v A ^ /10 2 GeV\ 1/2 

> 1.05 x — . (28) 



leV ~ V M *i 

This leads to a weaker bound than that obtained from the rare decays. The same constraints are valid for the 
case of ^2- Notice that one can satisfy these bounds easily. 

We emphasize that the most stringent constraints, which follow from K® — > imply that for 

MvpjUA >1600 GeV- eV it is possible to find values for the Yukawa couplings that are consistent with both 
neutrino oscillation data and the quantity Ak (hh). As we discuss below, one could expect to observe a 
significant number of LQ events at the LHC for LQ masses in the approximate range 400 < Mq, 1 < 1000 
GeV, or v/\ lying between about two and four eV. Since we have no independent handle on the value of v& 
nor any reason to preclude values lying in this range, we conclude that consistency of an LHC LQ discovery 
and the flavor changing decay — ► would imply a relatively large triplet vev in this model. 

2. fx — e Conversion Constraints 

Searches for charged lepton flavor violation (CLFV) in processes that conserve quark flavor can impose 
severe constraints on new physics. In the present case, searches for \i-e conversion in nuclei yield the most 
stringent constraints on the LQ couplings and masses. The tightest bound on the rate for this process has 



been obtained by the SINDRUM Collaboration||2l|] using the gold nucleus. The result for the conversion to 
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FIG. 5: Branching fractions for AT^ — > e + e~(black) and Jf£ — ► /i~/i + (blue) versus the lowest neutino mass for NH 
neglecting the phases, Mq, ± VA = 1200 (left), 1600 (right) GeV • eV. The horizontal lines represent the current 1 a 
experimental bounds Jl^l . 

capture ratio R^ e is 

A „ Tin- + A(Z,N) -> e~ + A(Z,N)] „ , „ 

Ra^e = Wr TTTi L 77^ \ \r Tvi < 7.0 X 10~ 13 90%C.L. . (29) 

^ e T\ji- + A(Z,N)-n/ + A(Z-l,N + l)] v ' 

This ratio can be expressed in terms of the reduced conversion and capture rates Rr^ e = ^conv/^capt 
where Q22D 

|pll|2|pl2|2 5 

^conv= 1 '' 4' " (VW + 2yW) 2 , (30) 

where V® and W n ) are overlap integrals involving the upper and lower components of the muon and 
electron wavef unctions and the proton (p) and neutron (n) number densities. Using the values for these 
integrals computed in Ref. [22], the experimental value for w ca p t and experimental limit on R^ e shown 
in Table. [Till we obtain 90% C.L. constraints on the LQ parameters as shown in Fig. [8] . There we give 
Rp!^, e corresponding to lightest neutrino mass in NH and IH when M^^a = 800 GeV • eV. We observe 
that in either case, there exist values of the LQ parameters for which the predicted ratio lies below the 90% 
C.L. limit (indicated by the horizontal line). Although the present limits exclude large portions of parameter 
space, they do not preclude the model entirely. Future searches for \i — > e conversion, such as the Mu2e 
experiment proposed for Fermilab or the PRIME experiment at JPARC, could improve the sensitivity be 
several orders of magnitude. A null result from these future experiments could rule out the LQ model 
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FIG. 6: Branching fractions for K® — > e + e~(black) and K\ — > /i~/^ + (blue) versus the lowest neutino mass for IH 
neglecting the phases, M^^a = 1200 (left), 1600 (right) GeV • eV. The horizontal lines represent the current 1 a 
experimental bounds ifl^i . 



Nucleus 






WcaptUOS" 1 ) 


^conv/ ^capt 


\fku 


0.0974 


0.146 


13.07 


< 7 x 10~ 13 



TABLE III: Values of the relevant parameters for fx — e conversion in Au 




under consideration here unless the product M^ 1 va is quite heavy. We will explore this possibility in a 
forthcoming study. 

Finally, we note that the leptoquark contribution to other lepton flavor violating processes, such as 
fi — ► ej and \i — > 3e, arise at one-loop level, and thus provide less stringent constraints. 



C. Other Constraints 

1. Oblique Parameters 

In this section we study the constraints implied by leptoquark loop contributions to electroweak precision 
observables (EWPO). Since the LQ-matter field couplings are tiny due to their dependence on Y u , we may 
safely neglect vertex, fermion propagator, and box graph loop corrections and concentrate on corrections to 
the gauge boson propagators. For sufficiently heavy LQs, one may characterize the leading effects of these 
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FIG. 7: Branching fractions for — > e /i + versus the lowest neutrino mass for NH (left) and IH (right) neglecting 
the phases, M^ua = 800 GeV • eV. The horizontal line represents the current 1 a experimental bound llql . 



corrections in terms of the oblique parameters, S, T, and U j23h . The U parameter is typically quite small 
and does not add a significant constraint in this case. We have computed the full set of LQ contributions to 
the gauge boson self-energy functions, and the resulting contributions to T and 5 are: 



a(M z )T 



1 



M? 



lW(0) - c 2 (fL zz (0) + yftz 7 (0) 



New 



N c 
8vr 2 i> 2 



2AM 2 



(3D 



4s|c| 



S = 



^Rei fl zz (M 2 z ) - Uzz(0) + 



New 
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where AM = M^ 2 — M& 1 , s (c) is the sine (cosine) of the weak mixing angle in the MS scheme, and "New" 
indicates the contribution from new physics (^1,2 )■ We have carried out renormalization at a scale \i = M z . 
The complete expressions for LQ contributions are given in Appendix |Bj while in Eqs. (13 1 132b we give 
approximate expressions in the limit that the mass splitting | AM| and Z-boson mass are small compared to 
the LQ mass M^. The latest global fit to EWPO yields for these parameters are S = -0.10±0.10(-0.08), 
T = -0.08 ± 0.11(+0.09), and U = 0.15 ± 0.11(+0.01), assuming a value for the SM Higgs boson mass 



Mh = 117 GeV (300 GeV) Biol . The most important constraint is from the T parameter which gives us a 



la bound of leptoquark mass splitting | AM| < 60 GeV. 
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FIG. 8: The rate of /i — e conversion in Au nuclei versus the lowest neutrino mass for NH(left) and IH(right) without 
any phases M^va = 800 GeV ■ cV. The horizontal line represents the current 90% C.L. experimental bound ifl^l . 

2. Collider Constraints and p-parameter 

The current constraint on the LQ mass, and M$ 2 , comes from the direct searches at the Teva- 



tron tl6J] 



Mp M > 250 GeV. (33) 

while for the values of the triplet vev one has, 

1 eV < v A < 1 GeV, (34) 

where the lower bound follows from the assumption that m u ~ 1 eV and that the Yukawa interactions are 
perturbative, while the upper bound is from the constraint on the electroweak p— parameter 

IV. LEPTOQUARK DECAYS AND NEUTRINO MASSES 



In this section we study the main features of the LQ decays taking into account the constraints on the 
neutrino masses and mixing. From Eq. (O one can conclude that the leptoquark (LQ) decays could be 
different in each spectrum for neutrinos and we will explore this interesting connection. 
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A. Main Features of the Leptoquarks Decays 

The decays of the leptoquarks into fermions — ► d{ and ^2 — ► di v, where d{ = d,s,b and 
ej = e, /i, r, are of most interest. The partial widths are given by 

1 1 1 2 1 1 ^ 3 1 2 

r(*j - d ie +) = Lj-Uf^, r(* 2 - dti/j) = y^M* 2 . (35) 

The widths are proportional to the Yukawa coupling squared, governed by vr? v jv\. 

The other competing decay modes come from the charged current gauge interaction. Depending on the 
mass splitting between the leptoquarks, one has the decays of the heavy leptoquark into the lighter one plus 
a real (W) or virtual (W*) gauge boson. For example, if AM = M^ 2 — Mq, 1 > 0, the ^ 2 decay rates for 
these processes are given by 

r(* a -iW) = ^Aia.r&.rU (36) 



327T7^ 



r(* 2 -+ v x w~* -+ *itt-) = 9 * v ^ AM3 fi Ji 



32xM*r V AM 2 ' 



(37) 



r(* 2 *ie-(M-)P.(P„)) = Jg^p, (38) 

r(^ 2 - - *i<rcO = 3T(^ 2 - * ie -(M-)P e (i> M ))- ( 39 ) 

where we have omitted the CKM factor associated with the Wqq' coupling in Eq. d39l ) since we are consid- 
ering the inclusive channels of qq', and where = Mj/M$ a . In Fig. [9] we show the branching fractions for 
the decays of ^2 versus its mass with M$ 1 = 250 GeV and va = 5 eV and 10 eV in (a) and (b), respec- 
tively, assuming that the Yukawa coupling is diagonal for simplicity. We see that once the mass difference 
is large and a real gauge boson channel is open, then it takes over the fermionic channels. In Fig.[l0]we plot 
the branching fractions versus va with = 250 GeV and Mp 2 = 300 GeV. Also in the second case: 
M^j = 400 GeV, and M^ 2 = 450 GeV. We see that for va < 10 eV, the Yukawa coupling is sufficiently 
large so that the fermionic channels become dominant. The features are the same for ^1 if it is heavier than 
*2- The presence of charged leptons in the final state (rather than missing energy associated with neutrinos) 
facilitates the use of collider observables to identify the connection between the leptoquark properties and 
the neutrino mass spectrum. Consequently, we will focus on the parameter region in which M^ 2 > M^ 1 . 
To further quantify this situation, we present the scatter plots in Fig. [TTJ in a plane of AM — v A under the 
condition T(^2 — ► <hvj) > Y(%>2 — > ^\W*). We see that if the mass difference is small enough, there is a 
broad range for the vev v A ~ 1 — 10 6 eV, although an eV value may be more natural for the light neutrino 
mass generation. 



4 AM 5 
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FIG. 9: Branching fractions of ^ decay for (a) va = 5 eV and (b) 10 eV, respectively, with = 250 GeV. 




FIG. 10: Branching fractions of * 2 decay for = 250 GeV, My 2 = 300 GeV and M-q, 1 = 400 GeV, M* 2 = 
450 GeV. 

B. Leptoquark Decays and Neutrino Spectra 

As discussed earlier, the lepton-fiavor contents of leptoquark decays will be different for each neutrino 
spectrum. Here, we study this issue in great detail. In Fig. [121 [13] and [14] we show the impact of the neutrino 
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FIG. 11: Allowed region of the parameter space when the leptonic decays of ^ are dominant for NH (left) and IH 
(right). AM versus va assuming = 250 GcV when 10~ 4 eV < mo < 10~ 2 cV, mo is the lightest neutrino 
mass. 

masses and mixing angles on the branching fractions of decaying into e, /x, r lepton respectively, with 
the left panels for the Normal Hierarchy (NH) and the right panels for the Inverted Hierarchy (IH). We first 
note that the absence of the decay channel be would indicate a NH with mi < eV, as seen in Fig. [12] 
In contrast, the be channel is the leading one in the IH. The s/i, dfi and dr, st channels have the leading 
branching fractions in the NH, but those channels alone do not seem to provide sufficient information to 
discriminate among the various possible neutrino mass spectra. If the lightest neutrino mass is above 0. 1 
eV, then this approaches the "quasi-degenerate" (QD) scenario of the neutrino mass pattern. The channels 
(be, s/i, dr) governed by the diagonal neutrino matrix elements reach to about 30% in this case. The 
interesting feature here is that d\i and the st channels vanish due to the unitarity cancelation of the mixing 
matrix. 

The predictions for the decays of ^2 taking into account the experimental constraints on neutrino mass 
and mixing angles are shown in Fig. [15] Because of the existence of the missing neutrino in the final state, 
we must sum over the contributing neutrinos incoherently. As one can see that, in the NH case the decay 
^2 — ► dv, sv is the dominant channel, and in the IH case ^ 2 — > b& is the leading one. In Figs. [16] and 
[IT] we plot the branching fractions of \&i the leading decay channels for NH and IH versus Mq, 1 without 
any constraints (red stars), and with Kl decay constraints and \i — e conversion constraints (green squares) 
when va = 4 eV. As noted earlier, present experimental results for these studies strongly constrain the LQ 
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FIG. 12: The branching fractions of — > d^e (d, = d, s, b) versus the lowest neutrino mass for NH (left) and IH 
(right) when all the phases vanish. Due to the symmetry as in Eq. (|20| |. the equal channels are also indicated in the 
parentheses. 

parameters, but do not preclude the possibility of scenarios in the LQ is sufficiently light to be discovered 
at the LHC. In the event of such a discovery, more precise measurements of B(Kl — > and more 

sensitive searches for [i — e conversion would provide a consistency test for this model. 

C. Impact of Majorana Phases in Leptoquark Decays 

Although the decays of the leptoquark are independent of the phases in the matrices K% and B, 
the unknown Majorana phases could modify the predictions for their decays. In this section we study the 
predictions for the LQ decays including the impact of the Majorana phases. It is important to note that the 
^2 decay branching ratios are independent of all unknown phases, including the phases in the K3 and B 
matrices and the Majorana phases, after summing over all final state neutrino flavors incoherently. This 
feature could allow us to probe the neutrino mass spectrum as well as the Majorana phases when combining 
the information from the \&i decays. To illustrate this point, we consider the limiting cases of mi ps 
(7713 ps 0) for the NH (IH), an approximation that we expect to hold when the mass of the lightest neutrino 
in either case is smaller than ~ 10~ 2 eV. For a heavier mass of mi 3 > 10 _1 eV, the situation approaches 
the quasi-degenerate, as we discuss below. 
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FIG. 13: The branching fractions of "fi — > di/j, (di — d, s, b) versus the lowest neutrino mass for NH (left) and IH 
(right) when all the phases vanish. Due to the symmetry as in Eq. (|20| |. the equal channels are also indicated in the 
parentheses. 

1. Normal Hierarchy with one massless neutrino (mi as 0) 

The leptoquark decay rates depend on only one Majorana phase $2 when mi w in the NH case. 
Taking S13 = for simplicity, one finds the expressions for the physical couplings: 



and 
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FIG. 14: The branching fractions of — > diT {d% = d, s, b) versus the lowest neutrino mass for NH (left) and IH 
(right) when all the phases vanish. Due to the symmetry as in Eq. (|20| |. the equal channels are also indicated in the 
parentheses. 

The behavior of the branching fractions for the dominant channels is shown in Fig. [18] We can see the rather 
weak dependence of the decay branching fractions on the phase $2 that can be also understood from the 
above analytical expressions. When the phase <I>2 = vr, one obtains the maximal suppression (enhancement) 
for the channels ^1 — ► sfi + and — ► dr + (^1 — ► dfi + , st + ), by a factor two at most. It is important to 
note that the braching ratio of the sum over all the quark contributions for each lepton flavor, s/x + + 
or dr + + st + , remains the same. Thus, from the observational point of view in collider experiments, the 
signals of /x+jet and r+jet will be unchanged, insensitive to the Majorana phase. 

2. Inverted Hierarchy with one massless neutrino (777,3 ~ 0) 

In the case of the Inverted Hierarchy with one massless neutrino one finds that all relevant decays depend 
on only one phase $1. The relevant couplings for S13 = are given by 

« /j^y e -^ 1(s 2 2+e -« lc 2 2)) (46) 
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FIG. 15: ^2 branching fractions versus the lowest neutrino mass for NH (left) and IH (right), independent of the 
phases. 
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FIG. 16: The branching fraction of leading channel 4*i — > s/i versus leptoquark mass for NH (left) and IH (right) 
without any constraints (red solid star), and with both Kl decay and fx — e conversion constraints (green empty 
square), va — 4 eV. 

In Fig. [19] we show the dependence of the branching fractions on this Majorana phase. The maximal 
suppression or enhancement takes place also when $i = ir. In this scenario the dominant channels also 
swap from one to the other. Unlike in the case of the NH where sfi <-> d/j, and st <-> dr at $2 = 7T, the 
interchanges of the channels occur for s/j,, dfx 6/x, and dr, sr <-> 6r, as well as be <-> cfe, se, when <I>i 
varies from to n. Once again, we notice that the branching ratio for the sum over all the quark contributions 
for each lepton flavor, be + + de + + se + or s^i + + d^ + + bfi + or dr + + sr + + 6r + remains the same. This 
fact makes the collider search for a leptonic channel rather phase-independent. However, if we require a 
6-flavor tagging in the event selection along with the lepton flavor identification, this could provide crucial 
information on the value of the Majorana phase $1. Specifically, we see the qualitative features in the IH 



BR(6e) < 



> BR(fyi), BR(6r) for $1 « 0; 
< BR(6/i), BR(6r) for $ x « vr. 



3. Quasi-Degenerate Hierarchy 

For a QD neutrino spectrum, the approximation of a nearly massless lightest neutrino no longer applies, 
and one must use the full expressions for the T 13 that are given in AppendixE] As one can readily appreciate 
from these expressions by setting m\ rs 7712 ~ m^, the branching ratios for the will in general depend on 
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FIG. 17: The branching fraction of leading channel "fi — > dr for NH (left) and — > be for IH (right) versus 
leptoquark mass without any constraints (red solid star), and with both K l decay and fj, — e conversion constraints 
(green empty square), va = 4 eV. 

both <I>i and <E>2- From a numerical scan over the neutrino masses and mixing angles, we find that the mild 
$2 -dependences of the branching ratios for <3?i = are quite similar to those of the NH, and that the strong 
<I>i -dependences of the branching ratios for $2 = are quite similar to those of the IH. Again, if 6-flavor 
tagging is effective along with the lepton flavor identification, one could hope to probe crucial information 
on the value of the Majorana phase $1. The only difference from the IH case is the numerical values for the 
branching ratio of the leading channels. 



D. Total Decay Width of Leptoquarks 

To complete our study about the LQ properties, in Fig. [20] we plot the total width (left axis) and decay 
length (right axis) for \&i versus v a for M^, 1 = 250 GeV in NH and IH. The total decay width is propor- 
tional to M$, 1 jv\. At the values va < 10 4 eV, its decay is prompt. This is the standard scenario for collider 
searches to be discussed in the next section. For larger values of va, the leptoquark can be long-lived in the 
detector's scale, making the collider signatures rather different. One may need to search for exotic hadrons 
or heavy charged tracks. Since these very large values of va are much higher than the light neutrino mass 
scale for our consideration of neutrino mass generation, we will not pursue this parameter region further. 
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FIG. 18: \&i branching fractions versus the Majorana phase $2 for the mi = scenario. 

V. SEARCH FOR LEPTOQUARKS AT THE LHC 

In this section we study the main production mechanisms of leptoquarks and their experimental signa- 
tures at the LHC. There exists extensive literature on this topic, treating both theoretical and phenomeno- 



logical considerations (25 



, 26, 27, 



s I29I. IboIi . ( 

goes beyond the existing studies in two respects: first we predict specific correlations between observ- 
ables involving different final state lepton flavors, and second we delineate the connections between these 
correlations and the light neutrino mixing angles and possible mass spectra. 

As is well known in the case of the leptoquarks, the leading production channel is via the QCD interac- 
tion 

q + q -» LQ + LQ (53) 
g + g LQ + LQ. (54) 

The pair production total cross section versus its mass at the LHC is plotted in Fig.|2T](the solid curve). The 
other unique channel is the single production via the Yukawa interaction 

g + q(q) LQ (LQ) + 1(1). (55) 



This cross section is rather small in our model due to the constraints from the neutrino masses, as we will 
discuss later. More importantly , the nature of the SU(2)l doublet $6 allows for associated production of 
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FIG. 19: branching fractions versus the Majorana phase $i for the 771,3 = scenario and $1 £ (0, 2ir). 
^1 and ^2 via W exchange, 

q(Pi) + q'(P2) - *i(fci) + n(k 2 ). (56) 

In terms of the polar angle variable y = f>\ ■ k\ in the partem cm. frame with energy yfs, the parton level 
cross section for the associated production is 

ay 32iV c sin 4 0v^ \ s ~ M w> 



where j3 = y[T — (mi + mj) 2 /s)(l — (m^ — rrij) 2 /s) is the speed factor of <I>i and ^2 in the cm. frame. 
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FIG. 20: Total decay width and length of versus va- 

The total cross section versus its mass at the LHC is plotted in Fig.[2l](the dotted curve), assuming = 
M^ 2 . As expected, the cross section for ^i^^ production is smaller than the QCD pair production 
or ^2^2 by about three orders of magnitude. 

We now turn to the signal observability at the LHC. In Sec. A and B, we are mainly concerned with 
the kinematical features for the signal and backgrounds. We will take the decay branching fractions of 
V&i, ^2 to be 100% to the corresponding channels under discussion. In Sec. C, we will devote ourself to 
the determination for the branching fractions. 



A. Pair Production of 



The QCD corrections to leptoquark pair production at the LHC have been studied B26I1 and a NLO K- 
f actor of order 1.5 — 1.9 for leptoquark mass range from 200 GeV to 1500 GeV is predicted. We apply 
the K-factor 1.5 to the processes in our numerical analysis. Assuming that V&i is the lighter leptoquark, its 
decay modes are 

*i -> <i;e+ (di = d,s,b; ej = e,fi,r). (58) 
We now explore the signal observability according to the different lepton flavors. 
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FIG. 21: Leptoquark pair production total cross section at the LHC versus leptoquark mass. The solid curve is for 
pp — > LQ + LQ via QCD, assuming /ip = fi R = My = M^ l . The dotted curve is forpp — > ^i 1 !^ + ^1 ^2 via 
exchange, assuming = AIy 2 = M*. 

We start from the cleanest channels with e, \i in the final state from \Pi decay. The signal consists of one 
pair of opposite-sign leptons of arbitrary e, \i flavor combinations plus two jets of d, s, b quarks. We employ 
the following basic acceptance cuts for the event selection J31I1 

pr(£) > 15 GeV, \q(£)\ < 2.5, (59) 
Pr(j) > 25 GeV, \r)(j)\ < 3.0, (60) 
A-Rjj, A%, AR U > 1.0 (61) 

where we require a large Ai? between the jet and lepton to resolve them because both ^is are produced 
nearly at rest implying that their two-body decay products are back-to-back. To simulate the detector effects 
on the energy-momentum measurements, we smear the electromagnetic energy, the electromagnetic energy 



and jet energy by a Gaussian distribution whose width is parameterized as 13111 

A.E a cat 



E ^/E/GeV 
Aff ahad 
E ~ y/E/GeV 



bcal, a cal = 10%, baa = 0.7%, (62) 
had, ahad = 50%, b had = 3%. (63) 
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The leading SM backgrounds to this channel are 

z*h*jj i + rjj, ti - £+rjj + Xt- (64) 

Although the background rates are very large to begin with, the signal and background kinematics are quite 
different. We outline the characteristics and propose some judicious cuts as follows. 

• For a few hundred GeV leptoquark decay, the leptons and jets in final states are very energetic. We 
tighten up the kinematical cuts 

pF*(l) > p? ax (i) > M 9l /A. (65) 

• To remove the jjZ* background, we veto the lepton pairs with opposite charges in the Z-mass 
window \Mg+£- — Mz\ > 15 GeV. This is a standard cut to remove the on-shell Z contribution. 

• To remove the ti background, we veto the events with large missing energy from W decay: jZ? < 
25 GeV. 

• In order to select the correct lepton and jet combination and reconstruct the leptoquark, we take 
advantage of the feature that the two leptoquarks have equal masses Mi 1 j 1 = Mg 2 j 2 . In practice, 
we take \Mi 1 j 1 — Mg 2 j 2 \ < M$ 1 /10. This helps for the background reduction, in particular for 

jn*/z*. 

The production cross section for the signal with the basic cuts (solid curve) and all of the cuts above 

(dotted curve) are plotted in Fig.|22] For comparison, the background processes of jZ* and ti are also 
included with the sequential cuts as indicated. Incidentally, the two background curves after the Z/W veto 
cuts coincide with each other. The backgrounds are suppressed substantially. There are also other sub- 
leading backgrounds, such as WZ, ZZ and W + W~ jj. Their production rates are much smaller compared 
to those mentioned above. For instance, including the decay branching fractions, we have 13211 

2 

a(WZ) ~ 17.25 pbx-x 6.72% = 0.77 pb, 
cr(ZZ) ~ 6.7 pb x 69.9% x 6.72% = 0.3 pb, 
a(W + W~ jj) ~ 0.2 pb with basic cuts. 

Therefore, these backgrounds become negligible after our selective cuts as discussed above. 

Finally, when we perform a signal significance analysis, we look for a resonance in the mass distribution 
of the Ij pair. The invariant mass of I + j is plotted in Fig. |23[a) for 400 GeV leptoquark production. The 
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FIG. 22: Production cross section of ^i^J with basic cuts and hard final states cut. Branching fractions for leptoquark 
decay are not included in this plot. For comparison, the background processes are also included with the sequential 
cuts as indicated. 



width of the distribution ~ ±20 GeV is governed by the detector resolution in our simulation. If we look at 
a mass window of \Mi 1 j 1: £ 2 j 2 — 1 < M^/20, the backgrounds will be at a negligible level. 

We would like to comment on the bb£ + £~ signal because the — > e + b channel play s an important role 
in distinguishing different neutrino mass spectra. The b-tagging rate is about 50% BIO . The signal is still 
sizable with two b jets tagged. 



2. T + i jj, T+T'jj 



Together with the e and fi produced in the decay, the r lepton final state from leptoquark decay can 
play an important role in distinguishing different neutrino mass patterns. Its identification and reconstruction 
are different from e, fi final states because a r decays promptly and there will always be missing neutrinos 
in t decay products. 

In order to reconstruct the events with rs we note that all the ts produced in the decay of few hundred 
GeV LQ are highly energetic. The missing momentum will be along the direction of the charged track. We 
thus assume the momentum of the missing neutrinos to be reconstructed by 



"^(invisible) = (track). 



(66) 
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FIG. 23: Reconstructed invariant mass of M(£j) for (a) £ + ( jj production and (b) i^r^jj production, with a 
leptoquark mass of 400 GeV. 



Identifying ^(invisible ) with the measured J^t, we thus obtain the r momentum by 

-Mi)- 



~Pt{t) = -pt(£) + Xt, P L (t) = p L (£) + 



The leptoquark pair kinematics are, thus, fully reconstructed. The reconstructed invariant masses of M(£j) 
and M(rj) are plotted in Fig. 1231 b). We see that M(rj) distribution is slightly broader as anticipated. The 
rather narrow mass peak of the £j system nevertheless serves as the most distinctive kinematical feature for 
the signal identification. 

For rrjj events with two r's, we generalize the momenta reconstruction to 



p (invisible) = «i p (tracki) + K2 ~p'(track2 



(67) 



The proportionality constants k\, K2 can be determined from the missing energy measurement as long as 
the two charge tracks are linearly independent. 

In practice when we wish to identify the events with r's, we require a minimal missing transverse energy 



> 20 GeV. 



(68) 



This will effectively separate them from the Ujj events. Another important difference between the leptons 
from the primary leptoquark decay and from the r decay is that latter is much softer. In Fig. [24] we show the 
Pt distribution of the softer lepton from the leptoquark and r decays in the events of Ujj, £rjj and rrjj, 
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FIG. 24: pt distribution of the softer lepton from the leptoquark and r decays in the events of ££jj, irjj and rrjj, 
for a leptoquark mass 400 GeV. 

B. Associated Production of ^i^^i 1 ^) 

The total cross section of ^1^2 + ^2 is shown in Fig. [2T]by the dotted curve. As discussed earlier, 
we focus on the fermionic decay ^2 — ► d*u. Correspondingly, we only consider the decays ^1 — > die and 
diii in order to effectively reconstruct the final state kinematics. The signal is thus 

#1*2, *i*2 -> i^j, vj (e = e,n), (69) 

consisting of one charged lepton and two jets plus missing energy. Once again, we consider the correspond- 
ing decay branching fractions to be 100%, and will leave the branching fraction determination to the next 
section. 

We employ the following basic acceptance cuts for the event selection 

pr(£) > 15 GeV, \r)(£)\ < 2.5, T > 40 GeV (70) 
Pr(j) > 25 GeV, \r](j)\ < 3.0, (71) 
ARjj, ARjt > 0.4. (72) 

To simulate the detector effects on the energy-momentum measurements, we use the same smearing param- 
eters as in the previous section. The irreducible SM backgrounds to this channel are 

W ± jj, W + W~, W ± Z -c e±jj +0 T . (73) 
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The Wjj background is by far the largest. To further optimize the signal observability, 

• we set additional cuts for hard lepton and jet 

pT x W > M*/4, pT x (j) > M*/4. (74) 

• The SM background events of £^J£t typically have the origin from W^. We thus veto the W 
boson with the help of a transverse mass cut Mt(£v) > 150 GeV. After this cut, the WW, WZ 
backgrounds are at a negligible level. 

• For the associated production of heavy particles like the two leptoquarks of several hundred GeV, 
the cluster mass of the system indicates the large threshold. We define 

M duster = ^M^ + C£p> T ) 2 + V(E^t) 2 +#t (75) 
and will impose a high mass cut to select the signal events M c i uster > 2M* . 

• The mass reconstruction for £j and vj can be very indicative. In order to select the correct lepton 
and jet combination to reconstruct fiwe define a transverse mass Mr by one jet and the missing 
transverse energy 

M T (ju) = yJ(Er(J) +Xt) 2 - (Mj) +Pt) 2 - (76) 

This variable has an upper bound Mjr, and is typically smaller than Mgj. We thus accept the momen- 
tum combination if M{ij\) > M T (j 2 v), M(£j 2 ) < Mt{j\v), or M(£j 2 ) > M T (Jiv), M(£ji) < 
Mt{j2v). However, there are two other possibilities that both transverse masses of the two combina- 
tions are smaller or larger than their invariant masses. In these cases, we select the combination with 
smaller \M(£j) — Mx{jv) \ as the right one. In Fig. [25] we show the reconstructed invariant mass and 
transverse mass distributions for a leptoquark mass 400 GeV for illustration. After selecting correct 
resonance, we impose transverse mass and invariant mass cuts for leptoquarks 

M T {jv) < Mvt, \M(£j) - AT* | < 20. (77) 

We list the signal and background cross sections with the consecutive cuts in Table. [TV] for the case of 
M<jr = 400 GeV. We see that the W jj background is at the O(0.1 fb) level. The production cross section 
of ^1^2 + ^1^2 versus its mass with basic cuts and all cuts are plotted in Fig. [26] We see that the cuts are 
highly efficient for the signal and the background processes of Wjj is under control up to M# ~ 1 TeV. 
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FIG. 25: Reconstructed invariant mass of M(£ ± j) and transverse mass Mt{]v) for a leptoquark mass 400 GeV. 

We would now like to consider the other potentially large backgrounds, namely, ti — > bb£ + £~ + J$t and 
leptoquark pair production via QCD ^i^J — > jj£ + t~. Demanding one isolated lepton with the basic cuts 
and vetoing the extra leptons in the range 

PT (£) > 10 GeV, \r](£)\ < 3.0, (78) 

the tt background is reduced by more than two orders of magnitude and the ifi^* background by about 
five orders of magnitude. The lepton vetoing is considered part of the basic cuts. With the additional Mt 
and M(£j) cuts, the backgrounds would be under control, as listed in Table. JV] 



C. Measuring Branching Fractions and Probing the Neutrino Mass Patterns 

So far, we have only studied the characteristic features of the signal and backgrounds for the leading 
channels and have taken the decay branching fractions of the leptoquarks to be 100%. As we presented 
earlier, the ^i, ^2 decay branching fractions and the light neutrino mass matrix are directly correlated. 
Measuring the BR's of different flavor combinations becomes crucial in understanding the neutrino mass 
pattern and thus the mass generation mechanism. In contrast, it has been customarily assumed 1 3, 3(1 that 
the leptoquark is flavor diagonal in each fermion generation. We list the leading reconstructible leptonic 
channels along with the predicted branching fractions in Table. |V] where the light quarks (d, s) can be 
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FIG. 26: Total cross section for ^1^2 + ^1^2 production after the basic cuts (solid curve) and all cuts (dotted curves) 
and the leading Wjj background after all cuts ("x" marks). 
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4.8 x 10 5 


2.3 x 10 4 
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9.5 


8.2 


0.42 


Ah = 600 GeV 


4.8 x 10 5 


7332 


6.7 


5.1 


4.7 


0.091 


Ah = 800 GeV 


4.8 x 10 5 


2945 


4.2 


2.0 


1.9 


0.018 


tt 














Ah = 400 GeV 


186 


11 


7.5 


0.6 


0.58 


















Ah = 400 GeV 


3.2 x 10~ 2 


2.4 x 10~ 2 


7.2 x 10~ 3 


7.2 x 10~ 3 


7.1 x 10~ 3 


5.4 x 10~ 3 



TABLE IV: ^1^2 + ^1^2 signal and irreducible background jjW and reducible backgrounds tt and \I>i\I>| for 
400 GeV leptoquarkmass. 



identified as jets, and the b quark may be flavor-tagged. We also associate these channels with predictions 
of the neutrino mass patterns. We reiterate some of the key observations in connecting the LHC signals of 
the leptoquarks to the neutrino parameters. 
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Signal channels 


Leading modes and BR 
Normal Hierarchy 


Leading modes and BR 
Inverted Hierarchy 


Leading modes and BR 
Quasi degenerate 


$! = $ 2 = 




e+be-b : (40 - 50%) 2 






fj+j w t+ j : 
any £ pairs: (40 - 60%) 2 


fi+j ps r+j : 
any £ pairs: (20 - 30%) 2 

fi + j e~b pa r + j e~6 : 
(20 - 30%) x (40 - 50%) 


e + b ps fi + s ps r + d : 
any £ pairs: (30%) 2 


$1 « 7T, $2 = 


independent of $1 


eb <-> ej, /ij <-» rj rb 
eb + ej, fij + /lb, rj + rb 
unchanged 


e& <-» ej, fis <-> /i&, rd <-> rb 
eb + ej, /ij + fib, rj + r& 
unchanged 


$1 = 0, $ 2 ~ T 


s/i «-> d/i, dr <-> sr 
jLtj, rj unchanged 


independent of <£> 2 


sfi <-> <i/i, dr <-> sr 
e6, /Ltj, rj unchanged 


*1*2(^1*2) 
$1 = $2 = 


(40 - 60%) x 100% 


e+bbv : (40 - 50%) x 50% 
fi+jbv: (20 - 30%) x 50% 


fi+sjD,e + bjv 2(30%) 2 
p+sW, e+bbv (30%) 2 


$1 PS 7T, $2 = 


independent of $1 


eb <-> ej, fij <-> /Lt6 


e£> <-> ej, /is <-> fib 


$1 = 0, $2 ~ 7T 


fi+j jD unchanged 


independent of $2 


eb, fis unchanged 



TABLE V: Leading fully reconstructible channels and the indicative ranges of their branching fractions. The light 
neutrino mass patterns of the NH and IH with the lightest neutrino mass mo < 10~ 2 cV are shown in the first two 
columns, including the effects of vanishing and large Majorana phases. The last column shows the Quasi degenerate 
neutrino mass pattern when mo > 0.1 cV. A jet of a light quark (d, s) is denoted by j. 

• ^2 decays are independent of any phases. They are thus robust to test the mass pattern. 

• ^1 decays to a lepton plus inclusive jets are also independent of any phases. They are robust to test 
the mass pattern as well. 

• These channels are not sensitive to Majorana phase <E>2- 

• The sensitivity to $1 can be significant if ^-flavor tagging can be exploited in the cases of the IH and 
QD. 

For illustration, consider first the cleanest channel, — ► fi + j fi~j (j = d, s). The number of events is 
written as 

N = L x a(pp *x*J) x BR 2 (^i -> fj+j), (79) 
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FIG. 27: Event contours in the BR— M*, plane at the LHC with an integrated luminosity 100 fb 1 for (a) ^i** — * 
fx~j, and (b) ^i^^ + ^1^2 — * H + j jv, including all the judicious cuts describe in the early sections. 

where L is the integrated luminosity. Given a sufficient number of events N, the mass of leptoquark is 
determined by the invariant mass of lepton and jet Mgj. We thus predict the corresponding production rate 
a(pp — > ^1^*) for this given mass. The only unknown in the Eq. ( p79T > is the decay branching fraction. 
We present the event contours in the BR— plane in Fig. |27[a) for 100 fb" 1 luminosity including all 
the judicious cuts described earlier, with which the backgrounds are insignificant. We see that the LHC 
has tremendous sensitivity to probe the channel ^1 — > /i+jet within and beyond the predicted branching 
fraction (35% - 60% in NH, 20% - 30% in IH, and 30% in QD) up to and beyond M 9l ~ 1 TeV. 

Since the ^2 decay is independent of the unknown Majorana phases, it is thus important to search for 
the ^^2(^1^2) signal. In Fig. 1271b). we show the event contours in the BR-M^ plane, for 100 fb" 1 
luminosity including all the judicious cuts described earlier. We see that with the estimated branching 
fraction for z/+jet, one can reach the coverage of about M$ 1 ~ 0.6 TeV or higher. 



D. Single Production of ^1 

The leptoquark pair production is well predicted by its SU(3) C and SU(2)l gauge interactions, and their 
decays studied in the previous sections are governed by the neutrino oscillation parameters, and all are 
model-independent. The drawback for this nice feature, however, is that it does not provide the model- 
dependent details, such as the Yukawa couplings or equivalently the vev v A . In order to further test the 
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FIG. 28: Total cross section for single production with 4 eV < va < 10 eV. 

theory, one needs to involve the Yukawa couplings in the production process. This is the single production 
of a leptoquark, as in Eq. (I55T ). 

The partem level cross section for this process is 

where t = (p^ — Pt^) 2 - Taking into account the neutrino oscillation constraints, in Fig. [28] we show the 
scatter plots for the total cross section of single production versus the lightest neutrino mass, for the 
NH case (left) and the IH case (right). The total cross sections for both cases are less than 1 fb, although it 
can become larger when reaching the quasi degenerate case for m v > 0.1 eV. To estimate the backgrounds, 
we start from the cleanest channels with e, /x in the final state from decay. The signal consists of one 
pair of opposite-sign same-flavor leptons and one light jet. We employ the same basic acceptance cuts and 
smearing parameters as the pair production. 

The irreducible SM backgrounds to this channel are 

Z*/j*j -»■ £ + rj, W + W~j -> £ + e~j +# T - (81) 

Once again, the signal and background kinematics are quite. We outline the characteristics and propose 
some judicious cuts as follows. 

• The leptons and jets from the leptoquark decay are very hard. We tighten up the px cuts 

pF*(£) > M*J4, P T x (j) > (82) 



s + t 



Ml 



+ 



(80) 



40 



<r(pb) 
cuts 


Basic 
Cuts 


cos 4>u 
> -0.5 


pfp , p J T CUt 

> 100, 100 GcV 


M z veto 
±15 GcV 


Mw veto 

Xt < 25 GeV 


M tj 
±20 GcV 


h*/z* 


309 


113 


9.3 


0.82 


0.82 


6.5 x 10~ 2 




0.87 


0.48 


5.8 x 10~ 2 


4.7 x 10~ 2 


5.6 x 10~ 3 


4.3 x 10~ 4 



TABLE VI: Backgrounds for a 400 GeV leptoquark single production. 



• Because the two leptons are largely back-toback for jZ* /j* background, we can set a cut for the 
angle between the two leptons in the transverse plane cos 4>u > —0.5. 

• To remove the Z background, we veto the lepton pairs with opposite charges in the Z-mass window 

\M t+t - -M z \> 15 GeV. 

• To remove the jWW background, we veto the missing energy from W decay and require 0? < 
25 GeV. 

• In order to select the correct lepton and jet combination and reconstruct the leptoquark, we make use 
of the feature that the lepton from leptoquark decay is harder than that produced with ^i, as seen in 
Fig.Hlta). 

• Finally, when we perform a signal significance analysis, we look for a resonance in the mass distri- 
bution of £j. The invariant mass of £+j is plotted in Fig. l29l b) with the event selection for 400 GeV 
leptoquark mass. The long tail in the high mass region is due to the wrong choice of the lepton. We 
propose to look at a mass window of \M(j — \ < Mq, 1 /20. 

We summarize the background rates with the above consecutive cuts, optimized for a Mq, 1 = 400 GeV 
signal in Table. |VT] Unfortunately, even after all these optimal cuts, the 7* /Zj background is still larger 
than the signal by more than an order of magnitude for both the NH and IH cases. Unless we are deep in 
the quasi degenerate situation when the Yukawa coupling is substantially larger, the single LQ production 
would not be observable above the huge SM background. We will not pursue this channel further. 

VI. SUMMARY AND CONCLUSIONS 

We have studied the properties of light scalar leptoquarks, = (vl/i , ^2) ~ (3, 2, 1/6), in the con- 
text of a simple grand unified theory where the neutrino masses are generated through the Type II seesaw 
mechanism. The SU(5) symmetry of this theory implies that the coupling of the leptoquarks to matter are 
governed by the neutrino mass matrix. We considered the low energy constraints on the model parameters 



41 




FIG. 29: (a) pr{&) distributions for the primary lepton and that from "Ji decay, and (b) invariant mass M(£j) for a 
400 GcV single leptoquark production. 

and explored the feasibility of testing the theory at the LHC with detailed parton4evel simulations for the 
signal and SM backgrounds. We focus mainly on the different scenarios where the semileptonic decays of 
leptoquarks are the dominant decay channels. We found the following interesting results: 

• Among the low-energy constraints coming from meson decays, lepton flavour violation, and elec- 
troweak precision data, the strongest bound is from the Kl — > /J + >~ decay, that leads to M^ 1 v A > 
1600 GeVeV. 

• The decays of the scalar leptoquarks are quite different in each spectrum for neutrino masses. Con- 
sider the most important decay. In the NH scenario, the leading decay channels are \i + j and 
r + j, while the e + j channel is absent. In the IH scenario, the e + b channel is the leading one, 
while fi + j and r + j are next. In the Quasi-degenerate case, the leading channels are shared by 

eb ps fj,s w rd. 

• The decays of ^ are independent of the unknown phases in the theory. One can use these decays to 
test the neutrino mass pattern without additional ambiguity. 

• The signal dependence on the Majorana phase $2 is very weak. However, the dependence on $1 in 
the case of IH and QD hierarchy can be strong, and determining the bi contributions could help probe 
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$i if we can effectively exploit 6-flavor tagging. Specifically, we have, for the IH or QD hierarchy, 



We reiterate that due to the symmetry in the neutrino mass matrix, the theory has the interesting 
prediction for the correlations BR(b/i) = BR(se) and BR(fe-r) = BR(de). We also note that the 
difference between the IH and QD cases is rather small, with lb branching fraction to be 30% in the 
QD case, compared to 40 — 50% in the IH case. 

• The leading production of the leptoquarks is via the QCD process, pp — ► ^i^^. We demonstrated 
that besides the clean di-lepton channels from e, fj,, the t final state can be effectively reconstructed as 
well. Even with only the clean channel /i + \i~ +jets, the signal observability can extend to ~ 1 
TeV and beyond. 

• Although the rate is smaller than for QCD pair production, the electroweak process of associated 
production ^i^^i^) can be very useful by making use of both ^i, ^2 decays simultaneously 
to identify the quantum numbers of the leptoquarks and to distinguish between the neutrino mass 
spectra. Even with only the clean channel of \x + .Xy+jets, the signal observability can reach about 

~ 0.6 TeV or higher. 

• We have found that the single leptoquark production via its Yukawa couplings is very small in both 
NH and IH spectrum. 

The discovery of these leptoquarks at the LHC could test this appealing scenario for neutrino mass 
generation at a fundamental level, and give us a hint about a possible candidate theory for grand unification. 
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APPENDIX A: EXPLICIT EXPRESSIONS FOR Tt AND T 2 



As seen from Eq. ©, we have 



and define 



and 



Ti = T 2 vl MNS Kl and T 2 = BK%Vp MNS 



diaq 

m v 
VA 



vtj|2 



X V 



A- 



/ e i/Jl \ / e iai N 



I? 



e i/32 
^e^ 3 J 



K 3 



e ia2 
e iaz 



22 2,222 2,222 2,22 2,22 2 
c 13 c 23 m 3 + c 12 c 23 s 13 m l + c 23 s 12 s 13 m 2 "+~ c 12 s 23 m 2 "+~ s 12 s 23 m l 



- 2ci2C 23 si2Si 3 S23 cos(<5)(m? - ml), 

v 2 _ 2 2 2,2 2 2,2 2 2 2,2 / 2 2,2 2 2\ 

r 2 — c 23 s 12 m l + c 13 s 23 m 3 + s 12 s 13 s 23 m 2 + c m c 23 m 2 + s 13 s 23 m lJ 

+ 2c 12 c 23 s 12 si 3 s 23 cos(5)(ml - m 2 2 ), 
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+ cj^e-^ms). (A12) 
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APPENDIX B: THE CONTRIBUTIONS OF LEPTOQUARKS TO THE OBLIQUE PARAMETERS 



We have computed the renormalized self-energy functions flyvil 2 ) m the MS scheme and find that 
ftw'(O) = for V, V = Z, 7. Thus, the T parameter is determined entirely by n^v^(O) an d reads 



-(M? + M? + — — — n In — | 

9 v 1 2 ; A/f 2 - A/f 2 A/f 2 



(Bl) 



where = My k (k = 1,2) and where the running parameters a and s on the RHS have been evaluated 
at the scale /j, = Mz- The S parameter, in contrast, requires inclusion of TLvv {M z ) 7^ for V, V = Z, 7, 
leading to 



= mcY i Q y \m\ 

+M 2 k In M 2 k - 2M 2 k F 1 (M|, M fe 2 , Ml) 



*{Mz)S = ^E[ M »^ M l)-T hMi 

n z k=l 



(B2) 



where the sum runs over the two LQ states having hypercharge Ylq and where 

F n (a,b,c) = / dx x n ln[o(l — x) + bx — cx(l — x) — ie] 
Jo 



(B3) 
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